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Abstract 

Some short history of few-body methods originated from the fa- 
mous Skornyakov-Ter-Martirosyan equation is given, including lat- 
est development of Faddeev formalism and Efimov states. The 3q 
system is shown to require an alternative, which is provided by the 
hyperspherical method (i^'-harmonics) which is highly successful for 
baryons. 



1 Introduction 

The Skornyakov and Ter-Martirosyan paper [Q] which appeared in 1956 marked 
the beginning of a new era in few-body physics, when a somewhat neglected 
part of nuclear physics was promoted to the successful domain of theoretical 
physics. As a result the few-body science has become a field accumulat- 
ing fast developing methods: L.D. Faddeev generalized the Skornyakov-Ter- 
Martirosyan Equation (STME) and has given a rigorous mathematical 
foundation for the theory of 3 particles [Q, many numerical methods have 
been introduced, for a review and references see [Q. As an immediate conse- 
quence of STME, a new effect was found in 1971, called the Efimov effect P], 
which is studied till now with respect to possible experimental consequences 

i- 

The STME and Faddeev technic is most useful when particles are nearly 
on-shell, so that e.g. 3-body results do not depend much on the potential 
shape, but rather described by the on-shell two-body t-matrix as it is for 
the quartet n — d scattering. The bound states of tritium and ^He provide 
another example, where the interaction at small distances (far off- shell) is 
important. To treat such systems an alternative method - the Hyperspherical 
Formalism (HF) (or iT-harmonics method) was developed and the system of 



1 



the Schroedinger-like equations was written |^. Its development was marked 
with many successful applications both in nuclear and atomic physical see e.g. 

1^. Recently it was understood that HF is probably the best suitable for 
systems with confinement such as 3 quarks, where the interaction is a three- 
body one, and confining so that the t-matrix formalism cannot be applied. 
Accuracy of HF as applied to the 3q system was found to be remarkably 
good 0, allowing for the 1% bias in the baryon mass [|12|. 

This talk is intended to demonstrate the physics of the 3-body system, 
and a qualitative analysis of two alternative approaches discussed above. 



2 The STME and Faddeev approach 

In the system of 3 equal-mass particles with arbitrary numeration one can 
introduce the total kinetic energy E and the momentum k in the pair (^,^ 

and the relative momentum p of particle 1, namely p = ^^g^-'' — |ki, k = 
k2-k3 

2 

The symmetric function of the ground state "^symm is expressed through 
partial w.v. 

'^symm = Pi) + ^/'(ksi, P2) + P3) (l) 

with the normalization condition 

\-^symm\^dVp = 1. (2) 

It is convenient to extract the free 3-body Green's function, introducing 

and the 3-body rescattering equation, equivalent to the summing the "bridge" 
Feynman diagrams (nonrelativistic) is 



X(k, p) = xo(k, p) - 2 y m p.2 + + p2 _ • 



(4) 



Here t{k,k',e) is the 2-body t-matrix, representing the "knot" in a bridge 
diagram, and 

^ ^(fc.? + Po.-g-ig)<^«(P+bo) 
Xok,p =-2m — — ^ , (5) 

P^ + Po + PPo 
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where is the 2-body bound state, while Po is the momentum of incident 
particle. 

Near the bound-state pole t-matrix can be written as 

fff :■;', + 0(r„) (6) 

[l-Kymla + ly Zme) 

where a = l/a,a is the scattering length and g{k,e) formfactor, g{0,0) = 1 
and g{k,e) fast decreases when k ~ I/tq and e ~ 

Let us assume now that the range of integration in (|^) is small p,p' <^ 
I/tq. Then one can insert in (^) with g = 1, and one gets- for the 3-body 
bound-state w.f. 



This is the STME for a 3-body bound state. The off-shell generalization 
of STME is the Faddeev equation(H). As it was correctly stated in the 
bound-state equation (0) cannot be used for tritium and ^He, since it has 
no lower bound for energy due to the Thomas theorem |T^. This can be 



easily understood rewriting (|^ in the form x = I Kxdp', and calculating 
the norm of K,\\ K ||= / (ip(ip(ii'(p, p'))^, which diverges logarithmically 
at large momenta, implying that there are formally infinitely many bound 
states. The physical situation corresponds to the cut-off form-factors g{k,e) 
present in K, which leads to the finite result for the norm || K \\. 

A specific situation occurs when the 2-body scattering length a is large, 
a ^ tq. Then the number of bound states lying between ( — ) and ( — ^ ) 
is approximately equal to 

iV ~ - In ^ (8) 

TT ro 

and when \a\ is increasing, \a\ oo, there appears an accumulation point of 
bound states (the Efimov effect 0). For 3 nucleons however < 1 and the 
effect is absent, but for three ^He atoms a = 104A, ro = 7 A and the effect is 
theoretically possible @]. 

Since the Efimov states are almost on-shell, it is convenient to calculate 
them using the 3-body unitarity and the N/ D method [|14|. Numerical results 



obtained (see Fig. 7 of |]T4[) support the estimate (|g) and yield the explicit 



position of levels near the energy threshold. 
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To conclude with the bound state equation (|^ it is interesting to study 
the properties of the bound wave function, e.g. the size of the bound system. 
Here one encounters an important difference between 2- and 3-body systems 
15| . Namely the 2-body loosely bound system with a small binding energy 
e, merQ <^ 1 has a radius of the order of r2 = -^^,^2 ^ r^. For the 3- 
body system the situation may be twofold. In case when a bound 2-body 
system exists as a subsystem and the 3-body bound state is close to the 
2-1-1 threshold, one has a quasi- two-body situation, whereas when 2-body 
bound subsystems are absent the size of the 3-body bound state is always rg 



however small binding energy is |T6[. Modern calculations of 3-body bound 
states in the framework of STME and its development - Faddeev equations 
are done for and ^He using systems of around 30 equations and exploiting 
realistic potentials describing NN scattering and bound states in the large 
energy interval (0-350 MeV). see e.g. the review in Unfortunately results 
of calculations yield significant underbinding of around 10-15%, may be due 
to three-body forces, which are not exactly known, hence final results are 



model dependent, see [0 for an example and references. 

We go now to the dN scattering which was also the topic of the primary 
paper The corresponding equations look like 

{y/3kyA-ME~at) /, 1 ^ _ -1 

^ a3/2(k,koj - ^2^^2^\^]^^^_ME~ 

_ r 47ra3/2(k .kp) dk' 

•' lu2 I u2 I L-lr' i\,i-c\(u2 u2\ 



(9) 



_ r 47ra3/2(K .Kp) dk . 

^(fc2+fc2+kk'-Af£;)(fc2-fc2) (27r)3' 

{Jzk^/i-ME-at) n 1 ^ _ 1/2 
PTfc^ a3/2(k,ko) - ,2+,2+kko-Afi? + 

. 47r{l/2ai/2 (k',ko)+3/2fei/2 (k'.kp)} ^k' 

^ (fc2+fc2+kk'-Af£;)(fc2^fcg) (27^)3 ' 

&l/2(k,ko) - fc2+;,2+kko-A/i? + 

, c 47r{3/2ai/2(k',ko)+l/2fei/2(k',ko)} ^k' 

^ (fc2+fc2+kk'-A/£;){fc2_fco-) (2'^)=* ■ 

Here 03/2 is the Nd quartet {S = 3/2) scattering amplitude, while 61/2 
and ai/2 are doublet {S = ^) amplitudes corresponding to the singlet and 
triplet last A^A^ interaction respectively. It is seen that the kernel for S = 
3/2 is mostly negative and allows for a faster convergence, in contrast to 
the doublet {S = 1/2) case. Numerical result for quartet scattering length 
(^3/2 = 5.1fm obtained in ||l[ is not far from experimental value [jl9l, whereas 
doublet scattering requires full off-shell calculation M. 
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3 Hyperspherical Method 



Heretofore the basic dynamics was assumed to be quasi-two-body (however 
the Faddeev technic allows for the full off-shell description) , in the sense that 
typical distance R between an interacting pair and a third spectator particle 
is large, R ^ Tq. However this situation is an exclusion, and not the rule, 
which can be understood from the representation of the w.f. through the 3 
body Green's function (^, r/ are Jacobi coordinates) 

V) = Mt ri) + J - 77 - r]')Vs(^', vM^'; r]')d^'drj'. (10) 

Here G{^, rj) = k = ^j2m\E\, p = ^/^T^, and V3 includes all 

interaction terms. The asymptotics of ip is given by G and is equal to 

^(^'*?)^~4' P^OO- (11) 

Hence the 3-body kynematics tends to concentrate all 3-body w.f. inside 
the interaction region of all 3 particles which generates small radius of w.f. 
even for barely bound 3-body states. (This is also true for -body systems 
> 3). In this situation any pair angular momentum l^j contributes to the 
total energy of the system an amount AE ~ ^'^2mJ^^^ which for the 3 nucleon 
system with tq ~ 1 fm and for l^j = 1 is of the order of AE ~ 50 MeV, while 
for the 3q system with m = ~ 0.3 GeV and Tq ~ 0.5 fm, AEq ~ 600 
MeV. 

Therefore it is advantageous to have a wave function with the minimal 
number of pair internal angular momenta for the given total momentum L. 
This basis is provided by hyperspherical functions {K - harmonics) due to 
the following properties: 

i) the solution of the condition lij"^ — 0,i ^ j — 1, ...N is given by the 
representation ^k=o — ^o(p), 

= ^ E (r. - r,)^ (12) 

i<j=l 

where all particles are assumed to have the same mass. 

ii) The function ^^(ri, ...r^) = UK{n)xK{p), where UK{n) = 

and Vk - harmonic polynomial contains excited angular momenta li, ...In-i 
the arithmetic sum of which is equal to K. 
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Therefore the basis corresponds to the minimal excitation of angular 
momenta and is advantageous for compact A^-body systems. Since as was 
explained the majority of such systems are compact, the Hyper spherical Ex- 
pansion Approach (HEA) formulated as a system of coupled integral or 
differential equations has proved to be very successful both for few-nucleon 
systems 0,[H1, where short-range correlations can be taken into account in 



the hyperspherical correlated basis (last ref. in IT^), and atomic physics 



It was understood afterwards [|TO|,[0 that the HEA works even better for 3g 
systems, since interaction there contains no repulsive core and confinement 
excludes two-body channels. 

Therefore already the lowest approximation with K = yields the 1% 
accuracy for the baryon energy []TU|, |TT], P^ . 

In this case the baryon state is characterized by the grand angular mo- 
mentum K and radial quantum number n = 0, 1, 2, which counts number of 



zeros of the w.f. in the p-space. A typical calculation was done in |T2|, and 
the result depends on only two input parameters: string tension a = 0.15 
GeV^ and = 0.4, while current masses of light quarks have been put to 
zero. The spin-averaged masses |(M„ -|- Ma) have been computed to elimi- 
nate effect of hyperfine splitting. 

To illustrate the simplicity of the method, let us quote the the equation 
for the for the dominant hyperspherical harmonics ipKip) = ^'^^ ; 



2fx dp 



+ Wkk{p)^k{p) = Ek^k{p) (13) 



where Wkk{p) is the sum of kinetic (angular) and potential energies, and 
/i is a constituent quark mass to be found below dynamically. The nonrela- 
tivistic appearance of this equation contains nevertheless the full relativistic 
dynamics, since p is the einbein field needed to get rid of square roots in the 
relativistic quark action. 

The explicit expression for Wkk is 

Wkk{p) = ^,+Vkk{p), d={K+^){K+l), (14) 

while Vkk{p) = {ux{^)yuK{^)), is the total potential V, including 2-body 
and 3-body parts, averaged over hyperspherical harmonics, which is done 
analytically. E.g. for the F-type 3q confining potential one has Vkk{p) = 



6 



1.58a p. It is remarkable that to find the eigenvalues with the 1% accu- 
racy one does not need to solve equation ([l3[), but instead is approximating 
Wkk{p) near the minimum point po by the oscillator well: 

WKKip)=WKKipo) + lip-pofW';,M, ^^|p=po=0. (15) 

The resulting eigenvalues are found immediately: 

EKn = WKK{po)+u^{n + ^), uj^ = W'^k/p. (16) 

The total baryon mass is calculated as Mku^p) = \p + Exnip), and finally 
p = Po is to be found from the stationary point condition ^^g^'*-^-* U=/^o — 0- 
This gives the constituent quark mass po = 0.957a/o" = 0.37 GeV and finally 
the baryon mass is MKnip^o)- The masses |(M„ + Ma) computed in this way 
are shown below in Table 1. 



Table 1 

Baryon masses (in GeV) averaged over hyperfine spin splitting for 
a = 0.15 GeV^ a, = 0.4, rrii = 0. 



State 


MKn+{^H)self 


{AH)coui 


nftot 


M*°*(exp) 


K = 0,n = 


1.36 


-0.274 


1.08 


1.08 


K = 0,n = l 


2.19 


-0.274 


1.91 


? 


K = 0,n = 2 


2.9 


-0.274 


2.62 


? 


K = L = l,n = 


1.85 


-0.217 


1.63 


1.6 


K = 2,n = 


2.23 


-0.186 


2.04 


? 



As it seen from the Table the calculated spin-averaged mass }^{M]^ + M/^) 
agrees well with the experimental average, the same is also true for lowest 
negative parity states with K = L = 1, which should be compared with 
1 \- states of N and A respectively. 

We also notice that breathing modes {n > 0) have excitation energy 
around 0.8 GeV while orbital excitations K = L = 1 have energy interval 
around 0.5 GeV. 

One of important advantages of HEA is that in the lowest approximation 
there is no need for numerical computations - as demonstrated above the 
result for the mass can be obtained analytically with 1% accuracy as can be 
checked by comparison with exact calculations, see [|l( 



12 . 
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To conclude, the on-shell approach of STME ( and its Faddeev general- 
ization) and the HEA are two alternatives which describe opposite physical 
situations. Their coexistence has played a very important stimulating role 
for the development of the few-body physics in the last four decades. 

The author is greatly indebted to L.N.Bogdanova for the help and useful 
discussions. 
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